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Strong confinement of charges in few electron systems such as in atoms, molecules and quantum 
dots leads to a spectrum of discrete energy levels that are often shared by several degenerate quantum 
states. Since the electronic structure is key to understanding their chemical properties, methods that 
probe these energy levels in situ are important. We show how electrostatic force detection using 
atomic force microscopy reveals the electronic structure of individual and coupled self-assembled 
quantum dots. An electron addition spectrum in the Coulomb blockade regime, resulting from a 
change in cantilever resonance frequency and dissipation during tunneling events, shows one by one 
electron charging of a dot. The spectra show clear level degeneracies in isolated quantum dots, 
supported by the first observation of predicted temperature-dependent shifts of Coulomb blockade 
peaks. Further, by scanning the surface we observe that several quantum dots may reside on what 
topologically appears to be just one. These images of grouped weakly and strongly coupled dots 
allow us to estimate their relative coupling strengths. 



The ability to confine single charges at discrete energy 
levels makes semiconductor quantum dots (QD) promis- 
ing candidates as a platform for quantum computation 
[U E] and single photon sources [3] . Tremendous progress 
has been made in not only understanding the proper- 
ties of single electrons in QDs but also in controlling 
their quantum states which is an essential prerequisite for 
quantum computation 0] . Single electron transport mea- 
surements have been the main experimental technique for 
investigating electron tunneling into quantum dots [S]. 
Charge sensing techniques using built-in charge sensors, 
such as quantum point contacts [Hj, complement trans- 
port measurements as lower electron tunneling rates can 
be monitored with even real-time detection being pos- 
sible j7]. It is instrumentally challenging to study self- 
assembled QDs via conventional transport and charge 
sensing methods due to the difficulty in attaching elec- 
trodes. Although progress is being made [51 fTUl [TTl [T2] 
these techniques have very small yield and therefore make 
it difficult to assess variation in QD electronic properties. 
Compared to typical QDs studied via transport measure- 
ments, in particular lithographically defined QDs, self- 
assembled QDs can be fabricated to have smaller sizes, 
stronger confinement potentials and a more scalable fab- 
rication process, all of which make them attractive for 
practical applications. 

In this paper, we focus on an alternative technique 
for studying QDs that is better suited for self-assembled 
QDs: charge sensing by atomic force microscopy (AFM). 
Charge sensing by AFM is a convenient method to study 
the electronic structure of QDs as nano-electrodes are not 
required and large numbers of QDs can be investigated 
in one experiment. Termed single-electron electrostatic 
force microscopy (e-EFM), this technique relies on the 



high force-sensitivity of AFM to detect the electrostatic 
force resulting from single electrons tunneling into and 
out of the QD. It was first demonstrated on QDs formed 
in carbon nanotubes [13l [14] , and later applied to self- 
assembled QDs [T5J [TB] and also to gold nanoparticles 
[TTJ [TS]. We focus on epitaxially grown self-assembled 
InAs/InP QDs in the few-electron regime. Using a dis- 
sipation model, we find compelling evidence for the ex- 
istence of electronic degeneracies (i.e. shell structure) 
by measuring an effective temperature-dependent level 
repulsion; although predicted for conductance measure- 
ments in 1991 ^19j, we believe this to be the first observa- 
tion of this effect. Further, we use the model to quantita- 
tively extract various properties of both individual QDs, 
such as the tunneling rates and charging energy and cou- 
pled QDs, such as the strength of coupling. 

We study uncapped self-assembled InAs QDs grown 
on a 20 nm InP tunnel barrier below which a two- 
dimensional electron gas (2DEG) is formed in an 
I1io.53Gao.47As quantum well. A dc-bias voltage, Vb, is 
applied to the 2DEG with respect to the grounded con- 
ductive AFM cantilever tip. Figure [T^ shows the sample 
structure and experimental setup. 

The AFM cantilever is driven at its mechanical reso- 
nance frequency, wo/27r ^^166 kHz, with constant oscil- 
lation amplitude [20]. The voltage drop, aVe (a < 1), 
across the tunnel barrier between the QD and the 2DEG 
is only a fraction of Vb, with a — a{x,y,z) being a 
function of the tip position. The tip-QD gap is tens 
of nanometers wide so that tip-QD tunneling is negligi- 
ble. We thus have a single-electron box setup: the elec- 
trochemical potential of the 2DEG, /i2DEGi with respect 
to the QD, //QD, is set by aVs and a negative bias in- 
creases the number of electrons, N, on the QD in integer 
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FIG. 1: e-EFM Ve-spectra and images, a, Schematic 
of the oscillating cantilever with pyramidal tip pushing elec- 
trons on and oft the QD when the mean bias voltage is just 
enough to lift the Coulomb blockade, b, y-Vs spectra taken 
at 4.5 K over upper QD shown in (e). Peaks in the spectra 
are always separated by the charging energy, but shells are 
additionally separated by AE. After passing each peak from 
right to left, the number of electrons, A'^, in the QD increases 
by one, with the N = 7 state after the leftmost peak. The 
energy difference between the first two peaks is 31 meV and 
the difference between peak 2 and 3 is 42 meV, so that if 2Ec 
between peak 2 and 3 is assumed to be 31 meV then AEsp 
= 11 meV. c. Topography of the InAs island with the ap- 
proximate locations of the QDs marked by X's. d and e. The 
simultaneously recorded frequency shift and dissipation im- 
ages of the structure in (c) at 4.5 K taken at —8 V. Scalebar 
is 20 nm. 



steps whenever the electrochemical potentials are aligned 
(called a charge degeneracy point). Tunneling between 
2DEG and QD is suppressed by the electrostatic energy 
cost, Eq, of adding or removing an electron to the QD 
except near these charge degeneracy points (Coulomb 
blockade) . The heart of the e-EFM technique lies in the 
fact that oscillations of the AFM cantilever modulate a, 
and hence are equivalent to an effective oscillating gate 
voltage applied to the QD. Thus, motion of the cantilever 
induces a modulation of N which will be slightly out-of- 
phase with the cantilever's motion (a result of the finite 
response time of electrons on the dot). The electrostatic 
coupling between QD and cantilever tip implies an elec- 
trostatic force proportional to N acting on the tip, the 
net result being both a frequency shift. Aw, and addi- 
tional dissipation, 7, of the cantilever [2T]. These effects 
are maximal at charge degeneracy points as here N can 
easily change in response to the effective oscillating gate 
voltage. 

Figure [Tja is an example of 7(Vb) at 4.5 K with the tip 
positioned over the upper QD imaged in Fig. [iji-e. The 
7(1^b) is equivalent to the energy addition spectra usually 
obtained from linear conductance or capacitance spec- 
troscopy measurements |22j . Similar to those measure- 



ments. Coulomb blockade peaks in 7 occur near charge 
degeneracy points of the QD. The peak positions are sug- 
gestive of the addition spectrum of a 2D circular QD with 
parabolic confinement potential; each peak is separated 
by twice the capacitive charging energy, 2Ec, with a fur- 
ther splitting between peaks 2 and 3 and between peaks 
6 and 7 corresponding to the energy difference between 
shells, AE. This type of shell structure has been previ- 
ously observed in InAs QDs [51 [51 [251 [M] . 

Figure [l]:-e shows the topography, Aoj and 7 images of 
an elongated InAs island. The peaks in the j-Vb spectra 
radially surround the QD center so that the ring furthest 
from the center corresponds to the first electron entering 
the QD; the rings themselves are contour lines of con- 
stant aVs- Multiple sets of concentric rings appearing 
in the Auj and 7 images indicate multiple QDs in the is- 
land. Such observations would not be as easily identified 
via other experimental means [23]. The tip-2DEG ca- 
pacitive force adds a large background in the Aw signal 
which locally varies due to topography [15l |26] , making 
it advantageous to focus on the 7 for image analysis. 

Figure [2^ shows 7 and Aui as a function of Vb with 
the tip positioned over the center of the upper set of 
concentric rings in Fig. [T^ at 30 K. While the energy 
addition spectra shown in Fig. [ij) and |2^ are consistent 
with the expected shell structure for a 2D circular QD 
with parabolic confinement potential, we obtain much 
stronger evidence of the shell structure from the tem- 
perature dependence of the peak positions. Theoreti- 
cally, a temperature-dependent shift of Coulomb block- 
ade peaks is expected whenever one has degenerate sin- 
gle particle levels, as predicted for the conductance peaks 
of a spin degenerate level [19]. The result is an effective 
temperature-dependent energy level repulsion: the peaks 
in each shell move apart as temperature is increased. Fur- 
thermore, our theoretical analysis suggests that this ef- 
fect should be enhanced in the tunneling-induced dissi- 
pation compared to the conductance due to an increased 
sensitivity to degeneracy. The same theory also leads to 
asymmetric lineshapes in the Coulomb blockade peaks, 
also enhanced compared to the conductance, but this is 
predicted to be small and could not be resolved in the 
present experiment. However, we were able to measure 
the temperature-dependent shifts of the peaks and find 
excellent agreement with theory. We are unaware of any 
experiments where these effects have been observed. 

We model the dissipation on the cantilever using lin- 
ear response and a master equation describing the charge 
state of the QD in the regime of weak coupling [THl [27] . 
Details of the approach are provided in the Supplemen- 
tary Information. Near a charge degeneracy point be- 
tween two nondegenerate single particle levels, the dissi- 
pation is [T51 [25] 
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FIG. 2: Temperature dependent shifts of Coulomb 
blockade peaks, a, The dissipation 7 (red) and frequency 
shift Awdip (blue) measured simultaneously versus bias volt- 
age, Vb, over the center of the upper set of concentric rings 
in Fig. [1^. Note that the parabolic background due to the 
capacitive force has been subtracted from Au to show AtJdip 
(see Methods). The lever arm, a, is determined by fitting 
each peak to equation ^ (black), b. The first two peaks 
from the right in 7 fitted to equation ([2| at different temper- 
atures (offset for clarity), c, The measured and theoretical 
separation between peaks 1 and 2 as a function of temper- 
ature where the sole fit parameter is the peak separation at 
zero temperature, 2i5ci =31 meV. A numerical calculation of 
this separation including the effects of the empty p shell is 
also shown (blue), d, The separation between peaks 3 and 6 
and peaks 4 and 5. Due to thermal broadening of the peaks, 
the positions of these peaks could only be determined up to 
30 K. The directions and magnitudes of the peak shifts as 
a function of temperature are indicated with arrows in (a), 
with larger arrows indicating greater shifts. 



where coq and fcg are the intrinsic cantilever oscillation 
frequency and spring constant, uj = luq + Alu is the mea- 
sured resonance frequency due to forces on the cantilever, 
r is the 2DEG-QD tunneling rate, and / = 1/(1 + 
exp{E/kBT)) is the Fermi function evaluated at i? = 
MQD - /i2DEG = ea{VB + Vq) (Vb = -Vb is the point of 
charge degeneracy). Lastly, A — —2Ec^ {I — a) ^g*'" 
is the sensitivity of the potential on the QD to the can- 
tilever motion and Ctip is the tip-QD capacitance. We 
stress that equation ([T]) applies to each dissipation peak 
independently: F and A are obtained separately for each 
peak from the data with no assumption of constant Eq. 

Equation ([T|) only takes into account single nondegen- 
erate levels. More generally, suppose we have a state with 

-I- 1 electrons, with risheii + 1 in the valence shell. If 
this shell has a degeneracy v, then near the charge de- 
generacy point between this state and the state with N 



QD electrons, the dissipation is: 

uj^A'^T (nsheii + l)(t^ - ?^shcn) /(I - /) 



i{Vb) 



where 



cj)= {y - nshc\l)f + ('^shcll + 1)(1 - /), 



(2) 



(3) 



and the tunneling rate F is assumed to be equal for each 
degenerate single particle level within a shell for a given 
peak [55]. Note that righeii is the number of electrons 
occupying the given shell and not the total number of 
electrons on the dot, N, and that because of the factor 
(j), 7(Vb) is no longer symmetric about its maximum. The 
different coefficients of / and 1 — / in reflect the asym- 
metry between electron addition and removal processes 
caused by degeneracy, and (f) — 1 corresponds to a nond- 



reduces to equation 
causes each peak in 



generate level for which equation ( 2 1 
([ij . This asymmetry in equation ( 2 1 
7(Vb) to be shifted in energy by an amount proportional 
to temperature. By fitting 7(Vb) (e.g. Fig. [2^) to equa- 
tion (|2| we extract a, allowing us to convert the Vb axis 
into energy. This is done for all of the thermally limited 
peaks, yielding a = 0.036±0.003. Figure[2]D shows 7(Vb) 
at different temperatures together with the fitted curves. 
The role of (j) is further elucidated by the relation 



(j)T — —2u!o 



Awdip 



(4) 



where Awdip is the size of the frequency shift dip due to 
the single-electron tunneling. The ratio in equation ^ 
defines an inverse timescale set by the relative in-phase 
and out-of-phase parts of the electrostatic force; this is 
simply F for a nondegenerate QD [T51 HH] , but modified 
by degeneracy through the factor (j>. Using equations 
([2| and ^ and the measured values of 7 and Awdip, we 
calculate the tunneling rates at the maxima of dissipation 
peaks 1-6, obtaining F/27r = 70, 90, 160, 180, 230 and 
330 kHz. As expected, F increases with increasing Vb 
as the height of the potential barrier between the 2DEG 
and the QD is reduced. 

After extracting the tunneling rates, we fit each dissi- 
pation peak using equation ^ and measure the spacing 
between peaks as functions of temperature from 4.5-30 K 
and from 78-95 K. We focus on the relative shifts between 
peaks as these are less sensitive to slight offsets in peak 
positions due to small changes in the tip-QD distance. 
The size and direction of each peak shift is different (see 
Fig.[2|L), in a manner that is completely captured by our 
model: the two peaks in the s shell shift apart, as do the 
four peaks in the p shell. The measured relative peak 
shifts of repelling pairs are shown in Fig. |2]:-d and com- 
pared to the theoretical shift from equation ^ with a 
single fit parameter, Eq- In addition, we expect that 
multiple shells, not just the valence shell, should play a 
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role at high temperatures where /cbT ^ Ai? is not satis- 
fied. For the relative shift of peaks 1 and 2 (Fig. ^) we 
plot a numerical calculation accounting for the possible 
occupation of the p shell, showing that the high temper- 
ature correction agrees well with the data. Figure 2d 
shows that the relative shifts between peaks 3 and 6 and 
between 4 and 5 are well described by equation ^ up 
to 30 K. Finally we note that there is an overall shift of 
the s and p shells toward each other that we believe is 
a consequence of strong repulsion of the p shell by the d 
shell, predicted to be 6-fold degenerate. 

We performed several checks to support our conclusion 
that the peak shifts result from degeneracy. First, the 
predicted shifts are unchanged if the relevant electronic 
levels are not perfectly degenerate, but rather split by an 
amount smaller than k^T; thus, the effect only requires 
approximate shell degeneracy as is the reality in imper- 
fect QDs. Second, if the level spUtting is larger than fceT 
we expect no temperature dependence of the peak shifts. 
Based on the observed peak shifts down to 4.5 K we 
thus conclude that the level splitting is smaller than this 
temperature, corresponding to roughly 0.4 meV. Lastly, 
the same degeneracy theory leads to small, but measur- 
able, shifts between the dissipation peak and the fre- 
quency peak corresponding to the same charge degen- 
eracy point. This is visible for the third peak (from 
the right) in Fig. |2]i in which the 7 and Au; peaks do 
not exactly line up, with the measured shifts compared 
to theory in Fig. S2. These dissipation- frequency shifts 
strongly support our model and rule out the alternative of 
a temperature-dependent renormalization of Eq or AE. 

Significant efforts are ongoing towards understanding 
and controlUng the properties of coupled QDs, in par- 
ticular double QDs or "artificial molecules" [SO]. The 7 
images that we obtain for double QDs are equivalent to 
stability diagrams which depict the charge state of the 
double QD system. This is because of the position de- 
pendence of the lever arm a(x, y, z) for each QD which 
results in two electrochemical potentials, /xqdi and ^qd2- 
Scanning the AFM tip at constant height and Vb corre- 
sponds to sweeping the /iQDi-MQD2 space through chang- 
ing tti and a2 even though only a single electrode (the 
tip) is being used. 

In a conventional stability diagram, lines of constant 
electrochemical potential for each QD are plotted as a 
function of two gate voltages. When the two QDs are 
coupled, intersection points are split into two points 
(triple points), showing avoided-crossings [30]. In the 7 
images the avoided-crossings are observed when the ring 
radii suddenly change at intersection points. Figure [3k,b 
show the same three QDs as in the lower part of FigT[l^, 
now imaged at —9 V and —7.6 V, respectively. Such 
avoided-crossings are highlighted in the circle and box 
in Fig. [3^, representing an example of weak and strong 
coupling respectively. 

We characterize the coupling strength by comparing 




FIG. 3: Imaging coupled QDs. a. Dissipation image show- 
ing the same 3 QDs as the lower half of Fig. ^ taken at a 
larger bias voltage, Vb ~ ^9 V. b. Dissipation image of the 
same region as (a) taken at Vb = —7.6 V. The 3 QDs are 
numbered in (a) but are easier to identify in (b) since each 
QD contains one electron. In (a) an example of weak coupling 
between QDl and QD2 is circled and an example of strong 
coupling between QD2 and QD3 is boxed. The possible mech- 
anisms are discussed in the text, c, Topography of two con- 
nected islands, d, Dissipation image taken at Vb = —8.0 V 
of the structure (c) . Each structure in (c) appears to have 
an associated QD. e, Dissipation image of region in (d) with 
many anti-crossings. Scalebar is 20 nm. 



the ratio of the change in ring radius to the separation 
between the first two rings (2£'c) for QD2. This method 
is only valid when both rings are far enough from the QD 
center that the voltage drop between them is approxi- 
mately linear. Following this procedure, the coupling of 
QD2 to QDl (circle) and QD3 (box) can be compared. 
While the change in radius of QD2 is approximately 0.10 
± 0.01 of 2Ec due to QDl, it is 0.46 ± 0.03 due to QD3, 
indicating a much stronger coupling between QD2 and 
QD3. We consider the former to be an example of weak 
coupling because the triple points are nearly joined. This 
is consistent with a small capacitive coupling between the 
two dots: the charging of one dot effectively gates the 
second dot, causing a sudden change in ring radius. 

Conversely, the boxed region in Fig. [3^ is an example of 
strong coupling since there is a large gap at intersections 
as in the triple points of a stability diagram. In Fig. |3]3, 
the same 3 QDs as in Fig. [3^ are imaged at smaller Vb . 
This image allows for a more intuitive explanation of the 
coupling. Consider the diagonal line from the center of 
QD3 outwards; initially, the AFM tip is over QD3 in the 
(^QD2, -^QD3)=(0,1) state, but takes a path into QD2 in 
the (1,0) state. The ability to go continuously between 
these states without going through (0,0) or (1,1) necessar- 
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FIG. 4: Imaging charge noise, a, Topography of QD. 
b, Dissipation image of (a) at 78 K. Middle sections of the 
outermost ring are missing, which indicates reconfiguration 
of nearby charge, c, Dissipation image of (a) at 4.5 K. The 
slight discontinuity between the upper and lower half of the 
rings indicate that the energy levels of the QD shifted to- 
wards ^2DEG in the lower half. Perhaps the cause of the 
"charge noise" is the emerging QD in the upper left corner 
of (c) which is in the vicinity of a small protrusion of the 
structure in (a), d, Dissipation image of a different QD at 
4.5 K with Vb = —8.0 V. During one horizontal scan line 
(slow scan direction is upwards), the confinement potential 
abruptly changes, e, Dissipation image of the same region as 
(d) at Vb ~ —9.0 V where coupled QDs are now present, f. 
Dissipation image taken at Vb ~ +6.8 V revealing no features 
as the QD is empty, g. Dissipation image taken immediately 
after (f) at Vb = —8.0 V. The original confinement potential 
was restored. Scalebar is 20 nm. (d)-(g) Share the same 
colorbar. 



ily indicates a large capacitive coupling between the dots. 
It also indicates evidence for an interesting charge trans- 
fer process as no dissipation is observed between circles. 
Lack of dissipation implies no change in the total dot 
charge; either there is a cotunneling process where two 
electrons simultaneously tunnel to and from the 2DEG, 
or there is coherent tunneling between the dots. 

Figure [3j;-e shows another example of coupled QDs at 
4.5 K. The InAs structure (Fig. 3:) contains coupled QDs 
as shown in the 7 image (Fig.Iji). Figure [3|3 zooms up 
on the region in Fig. [Sji showing many avoided crossings. 
Within this distance range from the QD centers, each a 
is approximately linearly dependent on the tip position 
so that scanning the tip more closely resembles sweeping 
two gate voltages, resulting in the image resembling a 
conventional stability diagram. Figure [3] also highlights 
how advantageous it is to have images in addition to the 
7-Vb spectra as the spectra alone will contain the peaks 
from nearby QDs which can be identified using the im- 
ages. 

Yet another advantage of our technique is the ability 
to spatially resolve the effects of background charge fluc- 
tuations; for self-assembled QDs, no other technique is 
capable of doing this. The AFM images show how the 



QD confinement potential is being infiuenced by charge 
reconfiguration. Figure |4] shows two such cases. The 
structure in Fig. |4^ shows a fiuctuation in electron pop- 
ulation due to nearby fiuctuations in the electrostatic 
background at both 78 K (Fig. |4|d) and 4.5 K (Fig. |4|:). 
The missing sections of the first ring in Fig. [4j3 indicate 
that the number of electrons in the QD is fluctuating by 
one in this region. Interestingly, depending on the scan 
direction (left to right or right to left) over the QD, the 
missing section may appear. Similar reconflguration was 
observed in the j-Vb spectra. While all the peaks ap- 
pear in the reverse Vb sweep (positive to negative), the 
flrst peak disappears in the forward sweep. Figure |4|l-g 
shows a more dramatic change. During the scan, a sud- 
den switch in the confinement potential occurs, leading 
to the transition from the single QD (Fig. |4|i) to a cou- 
pled double QD (Fig. [i]:). This double QD state could 



) by scanning 
4f). Although 



be switched back to the single state (Fig.Wa 
over the same area with a positive Vb (Fig. 
both of these changes are readily identified in the images, 
having a spectrum alone may cause confusion as is the 
case in conventional transport measurements. These ob- 
servations indicate that the QD confinement potentials 
are very sensitive to the electrostatic background and 
can be modified, or switched, controUably. 

Charge sensing with AFM can be used to investi- 
gate the electronic structure of single and coupled self- 
assembled QDs. The technique enables the quantitative 
extraction of the tunneling rate, charging energy, and 
the QD interaction energies; further, we have used it to 
perform the first measurement of temperature-dependent 
Coulomb blockade peak shifts confirming the shell de- 
generacy of the QD. The dissipation images proved es- 
pecially useful in analyzing multiple QDs and changes in 
QD confinement potential resulting from nearby charge 
fiuctuations. The images also revealed that what looked 
like a single QD structure topographically can actually 
contain multiple QDs. Additionally, the imaging capa- 
bility of AFM provides insight into the link between QD 
electronic structure and topography which is of great im- 
portance in developing electronic devices based on QDs. 



Methods 

The sample, grown by chemical beam epitaxy [5T| . 
consists of the following layers: 460 nm undoped InP 
grown on top of an insulating InP substrate, followed by 
a 10 nm Si-doped InP layer, 10 nm undoped layer, 10 nm 
Ino.53Gao.47As layer, 20 nm undoped layer and a 1.82 ML 
InAs layer that results in the formation of InAs QDs by 
Stranski-Krastanow growth. The QDs cover the surface 
with at a density of ~ 2.5 QDs per having diameters 
in the range of 30-95 nm and heights of 0.5-6 nm. The 
2DEG layer formed in the InGaAs well serves as a back 
electrode and an Ohmic contact to the 2DEG is made by 
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indium diffusion. 

Our home-built cryogenic AFM [32] includes an RF- 
modulated fiber optic interferometer |33| with 1550 nm 
wavelength for cantilever position detection. We coat Si 
AFM cantilevers (Nanosensors PPP-NCLR) with 10 nm 
titanium (adhesion layer) and 20 nm platinum. The can- 
tilevers typically have a 160 kHz resonance frequency and 
a quality factor between 100,000 and 200,000 at 4.5 K. 

All of the images were taken in frequency modulation 
mode [2D]. In this mode, the cantilever is self-oscillated 
at its resonance frequency with a constant amplitude. 
The frequency shift and dissipation were measured with 
a commercially available phase-locked loop frequency de- 
tector (Nanosurf, easyPLL plus). The topography im- 
ages were taken in constant frequency shift mode where 
a constant frequency shift is maintained by regulating 
the cantilever tip-sample distance using a feedback con- 
troller. The frequency shift and dissipation images were 
taken in constant-height mode with a typical tip height 
of 20 nm. Dissipation images are shown in Fig. SI as a 
function of Vb ■ More negative Vb results in adding more 
electrons to the quantum dot. Areas of increased dissi- 
pation mark 2DEG-QD tunneling events. Each time a 
ring is crossed when traveling towards the quantum dot 
center marks the addition of an electron to the dot. More 
details of the AFM images are listed in Supplementary 
Table 1. 

The amplitude of the cantilever excitation signal, Acxc, 
is provided as the dissipation signal from the Nanosurf 
oscillator controller. It is converted to units of 1/s via: 
10. ^cxc-A.xco ^ Acxco is the excitation amplitude indepen- 
dent of the tunneling process, in other words the back- 
ground dissipation. This conversion is independent of 
cantilever oscillation amplitude. Similarly, the signal is 
converted to units of eV/cycle by multiplying 

by the factor Eq 
lation amphtude [35 

The Aw— Vb spectra shown in Fig. [2^ was originally su- 
perposed onto a large parabolic background arising from 
the capacitive force between the 2DEG and cantilever 
tip. Over several volts, at typical cantilever-sample gaps 
of 20 nm, the curve can be fit with a single parabola. In 
Fig. |2ji this parabola was subtracted from the frequency 
shift data. 

The exact positions of the peaks (dips) in the dissipa- 
tion (frequency shift) are sensitive to the distance be- 
tween cantilever tip and quantum dot. In particular, 
slight changes in cantilever tip lateral position with re- 
spect to the quantum dot center can lead to slight shifts 
in the peaks as can be deduced from the images where 
the rings can have different spacing depending on loca- 
tion. The shift in peaks as a function of height, however, 
is linearly dependent over the distances used in this ex- 
periment (12-22 nm). We took the differences in peak 
positions in the data displayed in Fig. [2b to be caused 



"'^^ where a is the cantilever oscil- 

eQ 



by small height differences (sub 1 nm) and so the volt- 
age axis was rescaled to align the peaks with the data 
but the peak heights were not rescaled. The mean factor 
involved in the voltage rescaling is 1.011 with the most 
extreme factor being 1.088. The temperature data above 
22 K had thermally limited peaks for a cantilever oscil- 
lation amplitude of 0.4 nm, but needed to be reduced to 
0.2 nm at 4.5 K. The errorbars in Fig. ^ and d repre- 
sent how well the measurement over a single location was 
reproduced. 
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FIGURE DETAILS 



The tip height for the voltage spectra and all of the constant height images was 19 nm ± 1 nni, with the exception 
of Fig. 4d-g and Fig. SI where the height was ~ 23 nm. Additional image details are listed in Supplementary Table 
1. The acquisition time of the majority of spectra was 15 seconds. 
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FIG. SI: A series of constant height dissipation images, for the QD shown in Fig. 4d-g, for increasingly negative Vb- The base 
of the InAs structure is outlined with rectangular dashes and the highest area is outlined with rounded, more closely spaced, 
dashes. This QD is locaUzed near a high point in the structure, which is often observed. For increasingly negative Vb more 
rings emerge as the QD is populated with electrons. In these images the ring lineshape is broadened by the large cantilever 
oscillation amplitude of 0.33 nm at 4.5K with a tip-sample gap of roughly 23 nm. Note that the lateral position of the final 
image is slightly offset from the others as this image was taken at a later time in the experiment. Notice streaks appear in the 
same ring location indicating some nearby electrostatic influence. The scalebar is 20 nm. The same colorbar was used for each 
image, with all images but the last having a range of 0-0.85 Hz and the last 0-2 Hz. 



TABLE I: Experimental details of AFM images 



Fig. T (K) Auj/2n (Hz) Oscillation Amplitude (nm) Vb (V) Acquisition Time (min.) 
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DETAILS OF DISSIPATION WITH DEGENERATE SHELLS 



Here we outline the approach used to derive the general expression for the dissipation in Eq. (2). The charging 
Hamiltonian for small cantilever tip motion may be written 



He 



N 



{N-Afvf+ 1 + 



2DEG 



tip 



Hcfi + AHosc - J2 AnNz\N){N\, 



N 



|iV)(iV| 



(SI) 



is the dimensionless gate voltage, Ctip is the 



where \N) is a state with A'' electrons on the QD, Afy 
QD-tip capacitance and C2DEG is the QD-2DEG capacitance. In the second line, i?c,o is the cantilever-independent 
part of the charging Hamiltonian, AiJosc describes an electrostatic modification of the cantilever potential, and the 



dCtip 
dz 



We emphasize that Eq and A may be 



QD-cantilever coupling strength for given N is An = — 2i?c„-^ (1 — a) 
different for each electron added, as indicated by the index TV. From the second equality in equation (SI), N plays 
the role of a force on the cantilever. As a result, the dissipation and frequency shift may be found from the linear 
response coefficient AAr(w) describing the response of N to changes in z [27] . 

Consider the charge degeneracy point between N and N + 1 electrons on the QD. This may always be viewed 
as «shoii or nshcii + 1 electrons occupying a shell of degeneracy v (even for a nondegenerate single level, for which 
'^sheii = and — I). Neglecting interactions, the charge state with rtshcii ('^-shoU + 1) electrons in the shell is Z^n-fold 
(Dn+i-fold) degenerate, with 



"-shell 



D 



n+l 



"shell + 1/ ' 



(S2) 



where (') denotes a binomial coefficient. These arise simply from the different ways to put "shell or rishcii + 1 electrons 
into i> single particle states. Let Pnghcn.i be the probability to find nghcii electrons occupying the shell in configuration 
i, and Pusi-^au+i-j be the probability to find "shell + 1 electrons occupying the shell in configuration j. In general, these 
probabilities will satisfy the master equations |19j 



(S3) 
(S4) 



where Ti^j is the rate to add an electron to configuration i producing configuration j, and vice versa for Tj^i (note 
that these rates are nonzero only for configurations i and j that differ by the addition or removal of one electron). 
We calculate the rates using Fermi's golden rule. 

The master equations ( |S3| and ( |S4| may be solved for given values of and "shell, but the solutions are cumbersome 
for highly degenerate shells. To simplify the equations we assume that for a given charge degeneracy point (i.e. a 
single dissipation peak), the tunneling matrix elements from Fermi's golden rule are equal for all single particle states 
within the relevant shell. This is an approximation, since degenerate states may indeed have different wavefunctions 
leading to different tunneling rates. However, we expect the rates to be similar since the tunnel barrier between 
the QD and the 2DEG extends over the entire QD area, minimizing the effects of the spatial variations of different 
wavefunctions. Moreover, we checked that significantly unequal rates lead only to very small corrections in the peak 
shifts. For example, taking distinct rates for the two degenerate orbital states in the p shell, we find that rates 
differing by a factor of 2 lead to a correction of 1.5% for the shift of the 3rd dissipation peak (i.e. the 1st peak in 
the p shell). We thus neglect these possible differences here. Taking the rates to be equal we arrive at the simplified 
master equation for the total probability to find rishcii electrons in the shell, 



9t^«3hcll = ('^ - "shell) 



D 



(S5) 



where 



T+^Tf{E) , r_ = r [1 - /(£;)] 



(S6) 
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are the rates to add (+) or remove (— ) an electron to or from a single particle state, and / is the Fermi function. 
Note that the master equation for P„^,^^jj_|_i is not independent in our approximation of equal rates; this is a result of 



1. The stationary solution of equation (S5l is 

("shell 



P — 

^ "aholl 



1) 



Pn. 



hcll + 1 



{ly - nshcll) 



(1-/) 



(S7) 
(S8) 



where (j) is defined in equation (3). 

The quantity we need is the linear response coefficient AAr(a;). To find this, we assume that the cantilever is 
oscillating at frequency uj. This causes the chemical potential difference between the QD and the 2DEG to oscillate. 



E + 5e- 



and this leads to a change in the probabilities. 



^'"ahcll + l -Pnaholl + 1 + AAr(w)5e" 



Pn. 



(S9) 

(SIO) 
(Sll) 



Inserting equations (S7h(S11| into equation (S5l and linearizing in S, we solve for Xpf{u!). 
part yield the dissipative and conservative parts of the electrostatic force from (fco/'-^o)7 ~ 
{2ko/ujQ)Auj = A^^{X]^{uj)}. The dissipation for arbitrary degeneracy is given in equation (2). and for the frequency 
shift we obtain 



Its real and imaginary 
: -A^'i{XN{uj)}/uj and 



Auj = - 



2kQ k^T 



("shell + l)(i^ - y^shcu) 
Cj2 + (0r)2 



/(I-/)- 



(S12) 



Note that we recover the single level result [i.e. equation (1) for the dissipation] by taking v = \ and Tishcii = as 
expected. Finally, we point out that the temperature-dependent level repulsion discussed in the paper is contained 
in a symmetry of equations (2) and (S12l, from which we find that taking rishoii ^ v ~ "shell — 1 is equivalent to 
-E. 

The peak shifts of 7 and Aw are proportional to temperature and we can solve for the coefficients analytically. 
However, in general the coefficients are complicated and unenlightening. To show how the peak shifts depend on 
degeneracy, we provide the coefficients in the low and high frequency limits where they are greatly simplified. Note 
that our experiment is in the intermediate regime a; ^ F, so the peak shifts measured and calculated in the main text 
lie between these two limits. For 7, the peak shifts in the low and high frequency limits are 



AE. 



7, peak 



kBT 



'\n[d+ Vd(d+l) + l) 



as (w — > 0), 
as {ijj 00), 



where 



V — "shell 
"shell + 1 



1. 



(S13) 



(S14) 



For a nondegenerate level, d = Q and there is no peak shift at any frequency. For Aw the peak shifts in the same two 
limits are 



(S15) 



Comparing these limits, we see that the shell degeneracy results in a greater peak shift in 7 than in Aw. This is a 
direct consequence of equation (4), from which we see that, aside from an energy-independent prefactor. Aw differs 
from 7 by a factor of <f>. 

We measured the separation between the peak in 7 and the peak in Aw for each charge degeneracy point. This 




is shown in Fig. S2 for the third peak in Fig. 2a as a function of temperature and compared to theory with no fit 
parameters. As argued in the main text, this provides strong evidence that the observed peak shifts are indeed a 
result of shell degeneracy. 
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Temperature (K) 

FIG. S2: Difference between dissipation and frequency shift peak positions as a function of temperature for peak 3, compared 
to theoretical prediction with no fit parameters. 



